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Abstract 

By employing the QCD factorization approach, we calculated the next-to-leading order new 
physics contributions to the branching ratios, CP asymmetries, isospin and U-spin symmetry 
breaking of the exclusive decays B — ► Vj (V = K*,p), induced by the charged Higgs penguins in 
the general two-Higgs-doublet models. Within the considered parameter space, we found that (a) 
the new physics corrections to the observables are generally small in the model I and model III-A, 
moderate in model II, but large in model III-B; (b) from the well measured branching ratios and 
upper limits, a lower bound of Mh > 200 GeV in model II was obtained, while the allowed range 
of Mjj in model III-B is 226 < Mh < 293 GeV; these bounds are comparable with those from 
Qh! the inclusive B — » X s ^ decay; (c)the NLO Wilson coefficient Cfim^) in model III-B is positive 

and disfavered by the measured value of isospin symmetry breaking /S^{K*^) = (3.9 ± 4.8)%, 
but still can not be excluded if we take the large errors into account; (d) the CP asymmetry 



X 



Acp(B — * pj) in model III-B has an opposite sign with the one in the standard model (SM), 
which may be used as a good observable to distinguish the SM from model III-B; (e) the isospin 
symmetry breaking A(pj) is less than 10% in the region of 7 = [40 ~ 70]° preferred by the 
global fit result, but can be as large as 20 to 40% in the regions of 7 < 10° and 7 > 120°. 
The SM and model III-B predictions for A(pj) are opposite in sign for small or large values 
of the CKM angles; (f) the U-spin symmetry breaking AU(K*,p) in the SM and the general 
two-Higgs-doublet models is generally small in size: ~ 10~ 7 . 

PACS numbers: 13.20.He, 12.60.Fr, 14.40.Nd 
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I. INTRODUCTION 



As is well known, the inclusive radiative decays B — > X g 7 with q = (d, s) and the 
corresponding exclusive decays B — > V7 (V = K*,p) are very sensitive to the flavor 
structure of the standard model (SM) and to the new physics models beyond the SM and 
have been studied in great detail by many authors 0, 0, y|. 

For the inclusive B — > Xgj decay mode, the world average of the branching ratio |4j is 

B(B -> X sl ) = (3.34 ±0.38) x 1(T 4 (1) 

which agrees perfectly with the SM theoretical prediction at the next-to-leading order 
(NLO) El, |p , 5 1 B, ID] and puts perhaps most stringent bounds on many new physics 
models pjj, 11, U 13| where new particles such as the charged Higgs bosons may provide 
significant contributions through flavor changing loops. 

The exclusive decay B —> K*j has very clean experimental signal and low background, 
which was first observed by CLEO in 1992 jlij], and measured recently by BaBar and Belle 
with good precision 2, the world averages of the CP-averaged branching ratios are 
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B(B -> K 7) = (4.17 ± 0.23) x 10" 
B(B -> K*~j) = (4.18 ±0.32) x 10" 5 , (2) 

and have reached a statistical accuracy of better than 10%. The measurements of the 
Cabibbo suppressed B — > (p, 07)7 decays are difficult because the signal is about 20 times 
smaller, and the continuum background is about 3 times larger than the B — > if* 7 decay 
mode. Consequently, experiments have so far provided only upper bounds 3 HQ, but 



they will be surely measured at B factories in the near future. The currently available 
data as presented at the LP'2003 conference are summarized in Table H] 

When compared with the inclusive B — > X St cty decays, the corresponding exclusive 
B — > V'y decays are experimentally more tractable (specifically for B — > p7 mode) but 
theoretically less clean, since the bound state effects are essential and need to be de- 
scribed by some no-perturbative quantities like form factors and light-cone distribution 
amplitudes (LCDAs). 

In Refs.0, U^, the branching ratios and rate asymmetries of B — > decays 
were investigated in leading order (LO) and next-to-leading order (NLO) by employing 
the constituent quark model (CQM) jjl8]. In Ref.j2j]|, the exclusive B if* 7 decay 
was studied by using the perturbative QCD approach. Very recently, in the heavy quark 
limit nib 3> ^qcd, the decay amplitudes for the exclusive B — > (K*,p)j decay modes 
have been calculated in a model-independent way by using a QCD factorization approach 
"22! 23L i3 | , which is similar in spirit to the scheme developed earlier for the non-leptonic 



two-body decays of B meson |25|]. The NLO standard model predictions for the branching 
ratios, CP and isospin asymmetries, as well as the U-spin breaking effects for B —>■ K*j 
and B — » p7 decays are now available [13, H^, |2^, |23| . The new physics effects on 
isospin symmetry breaking and direct CP violation in B — > p7 decay have also been 
studied recently in supersymmetric models [28^ . 

In a previous paper, we calculated the NLO new physics contributions to the B° — B 
mixing and the inclusive B —>■ X s ^ decay from the charged Higgs loop diagrams in the 
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TABLE I: Experimental measurements of the CP-averaged branching ratios and/or CP violating 
asymmetries Aqp (at 90%C.L.) of the exclusive B — > V7 decays for V = K*,p and w. 



Channel 


CLE0[14] 


BaBar[15] 


Bellefl6] 


Average 


B(B -» K*° 7 )(10- 5 ) 


4.55 ±0.70 ±0.34 


4.23 ±0.40 ±0.22 


4.09 ±0.21 ±0.19 


4.17 ±0.23 


B(B -» K*+ 7 )(10- 5 ) 


3.76 ±0.86 ±0.28 


3.83 ±0.62 ±0.22 


4.40 ±0.33 ±0.24 


4.18 ±0.32 


B(B ^p°7)(10- 6 ) 


< 17 


< 1.2 


< 2.6 




B(B ^p + 7)(10" 6 ) 


< 13 


< 2.1 


< 2.7 




£>(-B — > tJ7)(10 _6 ) 




< 1.0 


< 4.4 




^ CP (i?^K* 7)(%) 


8 ± 13 ± 3 


-3.5 ±9.4 ±2.2 


-6.1 ± 5.9 ± 1.8 




Acp(B ^ K*+j)(%) 






±5.3 ±8.3 ±1.6 





third type of two-Higgs-doublet model (model III) and the conventional model II. In this 
paper, we calculate the new physics contributions to the branching ratios, CP asymmetries, 
and isospin and U-spin symmetry breaking of the exclusive radiative decays B — > (K*, p)j 
in the framework of the general two-Higgs-doublet models, including the conventional 
models I and II, and the model III. The QCD factorization method for exclusive B — > V'-f 
decays as presented in Refs. 22[ 23[ H| will be employed in our calculations. 

This paper is organized as follows. In Sec. II, we describe the basic structures of the 
general two-Higgs-doublet models, give a brief review about the calculation of B — > V'-f 
at NLO in QCD factorization approach in the SM and present the needed analytical 
formulas for the calculation of Wilson coefficients and physical observables. In Sec. Ill 
and IV, we calculate the NLO new physics contributions to the B — > K* / y and B —>■ p~f 
decay, respectively. The conclusions are included in the final section. 



II. THEORETICAL FRAMEWORK 



For the standard model part, we follow the procedure of Ref.[24] and use the formulas 
as presented in Refs. jiil. 2y]. The QCD factorization approach to the exclusive B — > 



V^y decays was applied independently in Refs. |22t 1231 |24J with some differences in the 
definition and explicit expressions of functions. We adopt the analytical formulas in the 
SM as presented in Refs. 0, HfjJ in this paper, since more details can be found there. 

In this section, we present the effective Hamiltonian and the relevant formulas for the 
exclusive decays B — > in the framework of the SM and the general two-Higgs-doublet 
models. 



A. Effective Hamiltonian for inclusive b — > 57 decay 



In the framework of the SM, if we only take into account operators up to dimension 6 
and put m s = 0, the effective Hamiltonian for b — > 57 transitions at the scale \i reads [24| 



n 



eff 



p=u,c 



C±Qi + C2Q2 ± ^ y CjQj 

j=3 



(3) 
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where A* = V* g V p b for q = (d,s) is the Cabibbo-Kobayashi-Maskawa (CKM) factor [29]. 
And the current- current, QCD penguin, electromagnetic and chromomagnetic dipole op- 
erators in the standard basis 1 are given by 2 



Qi = 


(sp) V - A (pb) v _ A , 


Ql = 




Qs = 


(sb) V - A ^2(qq) V -A , 


Qa = 


{s a bp)v-A ^2(q/3q a )v-A , 


Q 5 = 


(sb)v- A ^2(qq)v+A, 


Q 6 = 


{s a bp)v-A s ^2/{qpqa)v+A , 


Qi = 


g 

— —m h s a a^ v {\ + ^)b a F^ v , 

o7T 


Qs = 





where T a (a = 1, . . . , 8) stands for SU(3) C generators, a and (3 are color indices, e and 
g s are the electromagnetic and strong coupling constants, Q\ and Q 2 are current-current 
operators, Q 3 — Q G are the QCD penguin operators, Q 7 and Q 8 are the electromagnetic 
and chromomagnetic penguin operators. The effective Hamiltonian for b — > is obtained 
from Eqs.Q - © by the replacement s —> d. 

To calculate the exclusive B — > decays complete to next-to- leading order in QCD 
and to leading order in Aqcd/M b , only the NLO Wilson coefficient C 7 (/ib) and LO Wilson 
coefficients Ci(nb) with i = 1—6, 8 and /x& = C?(ra&) are needed. For the sake of the readers, 
we simply present these Wilson coefficients at the high matching scale My/ and the low 
energy scale fib = mb here, one can see [J, H| for more details. 

In literature, one usually uses certain linear combinations of the original Ci(fi), the so- 
called "effective coefficients" C cS (fi) introduced in Refs. 0, |3(} , in ones calculation. The 
corresponding transformations are of the form 

Cf(/i) = C,(Ai), (* = 1,...,6), (5) 

6 

Cf(iu) = CrM+^y&Oi), (6) 
i=i 

6 

Cf(iu) = C s (ii)+Y,ZiCi(f*), (7) 

i=i 

with y = (0, 0, 0, 0, -1/3, -1) and z = (0, 0, 0, 0, 1, 0) in the NDR scheme [jj, and y = 
(0, 0, -1/3, -4/9, -20/3, -80/9) and z = (0, 0, 1, 1/6, 20, -10/3) in the MS scheme with 
fully anticommuting 75 p. In order to simplify the notation we will also omit the label 
"eff" throughout this paper. 



1 There is another basis: the CMM basis, introduced by Chetyrkin, Mosiak, and Miinz 5] where the 
fully anticommuting 75 in dimensional regularization are employed. The corresponding operators and 

Wilson coefficients in the CMM basis are denoted as Pj and Zi in [24} . More details see 0, |3| . 

2 For the numbering of operators Q\ 2 , we use the same convention as Ref. 26] throughout this paper. 
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Within the SM and at the matching scale \i = My/, the leading order Wilson coefficients 

are 



C% M {M W ) = 1, (8) 
Cl SM {M w ) = 0,i = 2,..,6 (9) 



C° 7 , SM (M W ) = (10) 



Cl SM (M w ) = (11) 



with 



, , . 3x? - 2xi , -84 - 5a;? + 7x t ,„ _ 

AW = 4(^ ln " + • (12) 

^"'^ ^ffr (13) 

while the NLO results for C 7 (M W ) and C 8 (M W ) are 

! -16s«- 122s? + 80s? 1 fa* + 46s? - 28a? , 
<W Mwr ) = 9(i^l? + ln *' 



-1024 - 588^ - 2262:r t 3 + 32444 - 1364a; t + 208 
81(s t - l) 5 

16464 + 122054 - 107404 + 2509x t - 436 



\nx t 



486(x t - 1) 
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(14) 



-44 + 404 + 414 + x t 1 -174 - 314 2 

<W Mwr ) " 6(s* - l) 4 L ' 2(1 "4 )+ 2(s t -l)* lnXt 



-2104 + 10864 + 48934 + 28574 - 1994:r t + 208 

216(s t -l) 5 
7374 - 141024 - 282094 + 610s t - 508 
1296(:r t - l) 4 



In x t 



(15) 



where Xt = mf/rri^, and L%2(x) is the dilogarithm function. 

At the low energy scale ji = 0(m b ), the leading order Wilson coefficients are 



°i,sM = EM"'- for J = l,.-.,6, (16) 



i=i 



C? )S mM = ^Cl^Mw) + ^-^)Cl SM (M w ) + J2h l V a % (17) 

d i=i 

8 

c!,sM = cl SM (M w ) v ^ + M a * (18) 

1=1 
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in the standard basis, while 



Z lM = EVA for J = 1...6, (19) 



8=1 

= c° 7>s M, (20) 

Zjs-M = Cj SJ » (21) 

in the CMM basis. 

The NLO Wilson coefficient C-j(fib) at scale /i& = 0(mb) can be written as 

c 7 , s M = clM + iMc^M (22) 

with 



ClM = V%Cl SM (M w ) + ~ - V*) Cl, SM {M w ) 



/ 297664 is 7164416 M 256868 37 6698884 39 . n 

+ V 2S V 2S H V 23 V 2i Co q ,AM w ) 

V 14283 ' 357075 ' 14283 ' 357075 ' ' 



37208 
-I 1 n 23 

4761 1 ' 





C 7jSM (M w ) + £ (e^^(^) + /i + ft^rj* (23) 

i=l 



with 



. . xx 2 + llx-18 x 2 4x 2 - 16x + 15 ), 2, 2 

#(2: = — : — - + — ; M M --ni--, 24 

V ; 12(x-l) 3 6(x-l) 4 3 3' K J 

where r] = a s (Mjy)/a s (/z fe ), and the "magic numbers" a iy e^, /j, g iy Kji and hji, hi and hi 
are summarized in Table ITT1 

Using the central values of the input parameters as given in Table IHIl we find the the 
numerical results of the Wilson coefficients Cj(ra&) and Zj(m&) in the SM 

C°(m 6 ) = {1.1167, -0.2670, 0.0120, -0.0274, 0.0078, -0.0340, -0.3212, -0.1519} (25) 
Z°(m b ) = {-0.5339, 1.0280, -0.0055, -0.0727, 0.0005, 0.0012, -0.3212, -0.1519} (26) 

at the leading order, and 

C 7 {m b ) = -0.3212 + ,00112 = -0.3100 (27) 

C°(m b ) ACf LO 

at the next-to-leading order, the second term denotes the NLO QCD correction to Cj(rrib)- 



B. B — > Vj decay in the QCD factorization approach 

Based on the effective Hamiltonian for the quark level process b — > s(d)j, one can write 
down the amplitude for B — ► V'j and calculate the branching ratios and CP violating 
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TABLE II: The "magic numbers" appeared in the calculations of the the Wilson coefficients 
Ci{n) in the rare decay b — ► (77 with q = (d, s). 



i 


1 


2 


3 


4 


5 


6 


7 


8 


di 


14 

23 


16 
23 


6 
23 


12 

23 


0.4086 


-0.4230 


-0.8994 


0.1456 


k\i 








1 
_' 


1 

2 














k2i 








1 
2 


_ l 

2 














&3i 








1 

1 A 


1 




0.0510 


-0.1403 


-0.0113 


0.0054 










_ j_ 

1 A 
1*1 


_ 1 




0.0984 


0.1214 


0.0156 


0.0026 


km 














-0.0397 


0.0117 


-0.0025 


0.0304 


k&i 














0.0335 


0.0239 


-0.0462 


-0.0112 


hu 








1 


-1 






















2 
3 


1 

3 














h 3 i 








2 

63 


1 

27 


-0.0659 


0.0595 


-0.0218 


0.0335 


ha 








1 

21 


1 

9 


0.0237 


-0.0173 


-0.01336 


-0.0136 


hsi 








1 

126 


1 
108 


0.0094 


-0.01 


0.001 


-0.0017 


h§i 








1 

84 


1 

36 


0.0108 


0.0163 


0.0103 


0.0023 


ei 


4661194 
816831 


8516 
2217 








-1.9043 


-0.1008 


0.01216 


0.0183 


fi 


-17.3023 


8.5027 


4.5508 


0.7519 


2.0040 


0.7476 


-0.5358 


0.0914 


9i 


14.8088 


-10.809 


-0.8740 0.4218 


-2.9347 


0.3971 


0.1600 


0.0225 


hi 


2.2996 


-1.0880 


3 
7 


1 

14 


-0.6494 


-0.0380 


-0.0185 


-0.0057 


hi 


0.8623 











-0.9135 


0.0873 


-0.0571 


0.0209 



TABLE III: Values of the input parameters used in the numerical calculations |3ll. EES- For 
the value of Fjc*, we use the lattice QCD determination of Fk* = 0.25 ± 0.06 33] instead of the 
result Fk* = 0.38 ±0.06 as given in Ref-j^]. The smaller value of Fjc* gives a better agreement 
between the SM predictions and the data. Rb = \J p 2 + fj 2 , and A, A, p and fj are the ordinary 
Wolfenstein parameters of the CKM mixing matrix. 



A 


A 




7 


G F 




0.854 


0.2196 


0.39 ± 0.08 (60 ± 20)° 1.1664 x lO^GeV' 2 


1/137.036 


a s {M z ) 


mw 


m t 


A^L 

MS 


m c (m b ) 


m u 


0.119 


80.42 GeV 


174.3 GeV 


225 MeV 


1.3 ±0.2 GeV 


4.2 MeV 


Jb 




m Bd 


m b (m b ) 


1~B+ 




200 MeV 


(350 ± 150) MeV 5.279 GeV 


4.2 GeV 


1.671ps 


1.537pa 


F K * 


Ik* 


fk- 


m K * 


af 


af 


0.25 ±0.06 


230 MeV 


185 MeV 


894 MeV 


0.2 


0.04 


Fp 


fp 


J P 


m p 


a i 


«2 


0.29 ±0.04 


200 MeV 


160 MeV 


770 MeV 





0.2 
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asymmetries once a method is derived for computing the hadronic matrix elements. One 
typical numerical result obtained by employing the constituent quark model is 

B(B -> K*j) « 5 x 10" 5 (28) 

at both LO and NLO level |2(|. Although this theoretical prediction is in good agreement 
with the data numerically, but the hadronic models used in 18, 0, 2(| did not allow a 



clear separation of short- and long-distance dynamics and a clean distinction of model- 
dependent and model-independent features. By using the QCD factorization approach [i^, 
l23l . |24( , one can separate systematically perturbatively calculable hard scattering kernels 
( 3/ and 3/ ) from nonperturbative form factors and universal light-cone distribution 
amplitudes of B, K* and p mesons. The higher order QCD corrections can therefore be 
taken into account consistently. 

In this paper, we calculate the new physics contributions to the exclusive decays 
B — > K*^j and B — ► /ry in the general two-Higgs-doublet models by employing the 
QCD factorization approach. We will always consider the decay widths or branching 
ratios averaged over the charge conjugated modes with an obvious exception of the CP 
asymmetries. 

In QCD factorization approach, the hadronic matrix elements of the operators Qi with 
i — 1, ... ,8 for B — > V'-f decays can be written as [2J] 



{Vl{e)\Qi\B) 



(0)3/+ f didvTl\^v)^ B {i)^ v {v) 
Jo 



■ e (29) 



where e is the photon polarization 4-vector, F B ^ V is the form factor describing B — ► V 
decays, $s and <3V are the universal and nonperturbative light-cone distribution ampli- 
tudes for B and V meson respectively 3 , v ( v = 1 — v ) is the momentum fraction of a 
quark (anti-quark) inside a light meson: l{ = vk + and l~£ = vk + while = (k + , k~, k±) 
is a four vector in the light-cone coordinator, £ describes the momentum fraction of the 
light spectator quark inside a B meson: l + = with £ = 0(AQ C D/ m b), and 3/ and 3/ 7 
denote the perturbative short-distance interactions. The QCD factorization formula (J29|) 
holds up to corrections of relative order KqqdI m b- 

In the heavy quark limit, the contributions to the exclusive B — > V'-f decay can be 
classified into three classes 4 : 

1. The "Type-I" or "hard vertex" contributions include (a) the contribution of the 
magnetic penguin operator Q 7 described by form factor F B ^ V , which is the only 
contribution to the amplitude of B — > at the LO approximation; and (b) the 
0(a s ) contribution to the hard-scattering kernels 3/ from four-quark operators Qi...e 
and the chromo magnetic penguin operator Q s . 

2. The "Type-II" or "hard spectator" contributions include the 0(a s ) contribution to 
the hard-scattering kernels 3/ from four-quark operators Q1...6 and the chromo- 
magnetic penguin operator Q 8 . 



3 For explicit expressions and more details about $b and $y, one can see Refs.|22l l25j and references 
therein. 

4 For more details of various contributions and the corresponding Feynman loops, see for example Ref.|26| 
and references therein. 
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3. The "Weak annihilation" contribution, which is suppressed by one power AQCD/^b 
when compared with the Type-I and II contributions, and the dominant annihilation 
amplitudes can be computed within QCD factorization. 



Combining all parts together, the decay amplitude to 0(a s 
decay takes the form of 

A (B^ y 7 ) = ^£ Rv ( Vl \Q 7 \B), 



for exclusive B 



with 



R v = A? K(V T ) + a% nn (Vy)] + A<«> [o?(V 7 ) + a c ann (V 7 )) 



(30) 



(31) 



where q = s for V = K*, q = d for V = p, and a% (p = u, c) denote the hard vertex and 
hard spectator NLO contributions 



a^Vj) = CM + 



a s {ii)Ci 



An 
An 



zK»)Gi(z p ) + z °M G i 

i=3...6,8 



i=l,2 



J=3...6,8 



(32) 



where z„ 



m 2 /m 



a/0.5/1, CV = 4/3, the Wilson coefficients can be found in 
previous subsection, the explicit expressions of the functions Gi and Hj can be found in 
Ref . [jfil] and in Appendix |XJ The functions a^ nn and a^ n in above equation denote the 
weak annihilation contributions and take the form of 26] 



ai n (;r 7 ) = 

< nn (X*° 7 ) = 

< nn (ir- 7 ) = 

<nn{K*-l) = 

for B — > ff* 7 decays, and 

«an„(P°7) 
«ann(P"7) 
<nn(P°7) 



Qd [a A b K * +a 6 (df + df )], 

Q« [ai&^* + a 4 b K * + a 6 (-2df * + df)] 
Q u [a 4 b K * + a e (-2df + df )] , 



[-a 2 & p + a A b p + a 6 (d£ + dg)] 
Q u [a^ + a 4 fe p + a 6 (-2d£ + 
Qd [a 4 & p + a 6 (d£ + d£)] , 
Q u [a 4 &" + a 6 (-2d£ - 



(33) 



(34) 



for 5 — > p 7 decays, where Q u = 2/3 and Qd = —1/3 are the electric charge of up and 
down quarks, <2j denote the combinations of LO Wilson coefficients 



Ol,2 

o. i 

«6 



i /^O 

°i,2 + 3 2 4' 

°6 + q°5" 



(35) 



9 



And finally the functions b v , d\ and d\ are [2 

b v = 2tt 2 f B m v f v 
F v m B m b X B '' 



(36) 



d ^ = - Aj f¥t^ TaX+a " + -^- (37) 

v ry mBrrib 



The values of all parameters appeared in above two equations can be found in Table IIII1 
One special feature of the B — > /ry decay is that the weak annihilation can proceed 

through the current-current operator with large Wilson coefficient C\. Although the 

annihilation contribution is power-suppressed in l/m&, but it is compensated by the large 

Wilson coefficient and the occurrence of annihilation at tree level. 

From the decay amplitude in Eq.()30jh it is straightforward to write down the branching 

ratio for B — > V7 decay 

G \am\ml ( _ 3 , D |2 ^ 2 
32tt 4 I m? 



B(B — > Vj) =tb ~" a I 1 — -y |fl v r^|Fv| a , (3 



where function i?y has been given in Eq. (j3~T]) . and cy = 1 for V = K*, p and cy = l/y/2 
for V = p°. The branching ratios for the CP-conjugated B — * V"7 decay are obtained by 
the replacement of X p g ^ — > Xp in function i?y. 



C. Outline of the general 2HDM's 

The simplest extension of the SM is the so-called two-Higgs-doublet models In 



such models, the tree level flavor changing neutral currents are absent if one introduces an 
ad hoc discrete symmetry to constrain the 2HDM scalar potential and Yukawa Lagrangian. 
Let us consider a Yukawa Lagrangian of the form|34| 

U = n'ijQi.i.oi'',!.- + n', , ! Q,.iP\i\,n + £%Qi,LfaU j)R + ^Q i>L <p2D j)R + H.c, (39) 

where 0j (i = 1, 2) are the two Higgs doublets, 4>i,2 = i^^i 21 Q%l {U^r) with i = (1, 2, 3) 
are the left-handed isodoublet quarks (right-handed up-type quarks), Dj R are the right- 
handed isosinglet down-type quarks, while rf^ and £ 4 '• = 1,2,3 are family index ) 
are generally the nondiagonal matrices of the Yukawa coupling. By imposing the discrete 
symmetry 

0i -> -0i, 02 -> 02, A -> -A, ^ -> T^< (40) 

one obtains the so called Model I and Model II. 

In model III j^HHH, the third type of two-Higgs-doublet models, no discrete symmetry 
is imposed and both up- and down- type quarks then may have diagonal and/or flavor 
changing couplings with 0i and 02- As described in (24), one can choose a suitable basis 
to express two Higgs doublet 0i and $2 and define the mass eig enstates {H ± , H , h°, A 



After the rotation of quark fields, the Yukawa Lagrangian of quarks are of the form 

Zy 11 = vuQi,Lfau jtR + ///;g,,o, /;,/,> + BQ^fcu^R + Bq^d^ + h.c, (4i) 
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where r\(- D correspond to the diagonal mass matrices of up- and down-type quarks, while 
the neutral and charged flavor changing couplings will be |34| 

in'eutral = ') Charged = ^^CKM, £diarged = VcKm(, D , (42) 

with 

rj n q ./frurfu 
& D = *\ % , (43) 
3 V2M W 1 



where Vckm is the CKM mixing matrix [23], i,j = (1, 2, 3) are the generation index. The 
coupling constants Ay are free parameters to be determined by experiments, and they 
may also be complex. 

The two-Higgs doublet models have been studied extensively in literature at LO and 
NLO level 0,13,0123,121123,121123,53,51 53,53 and tested experimentally gj. 



For the model I, the new physics corrections to physical observables are usually very 
small and less interesting phenomenologically. The model II, however, has been very 
popular, since it is the building block of the minimal supersymmetric standard model and 
may provide large contributions to the mixing and decay processes of K and B meson 
systems. The most stringent constraint on model II may come from the inclusive B —>■ X s j 
decay. From the experimental measurements and currently available studies at NLO level 
' SHU! Ill, one get to know the following main features of the conventional models 



I and II, and the model III 

• For the model I, no bound on M# can be obtained from B — > X s ^ j44j ]. since 
the charged Higgs loops interfere destructively with the SM penguin diagrams and 
decouple for large tan (3. 

• In model II, the charged Higgs penguins interfere constructively with their SM 
counterparts, and thus always enhance the branching ratio B(B —>■ X s j). The 
measured mass splitting AM Bd = 0.502 ± 0.007ps _1 and the decay rate B(B — > 
Xg^) = (3.34 ± 0.38) x 10~ 4 leads to strong bounds on both the tan/3 = v 2 /vi and 
the mass M# jl3l . 44j ]. The typical bounds at NLO level as given for example in 



Ref. jil 



are 

tan/5>0.6 (44) 

and 



M H > 300GW (45) 

for any value of tan/5, and the tan/3 dependence of the lower bound saturates 
for tan/3 > 5. This NLO lower bound on M# is much stronger than the direct 
experimental bound M# > 78.6 GeV |3lJ and the bound from other observables, 
such as Rb and B — > r decays [8(- 
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For the model III 5 , the charged-Higgs loop diagrams can provide significant contri- 
butions to B° - B mixing, the inclusive B — > X s ^ decay and many other physical 
observables 34l 41 . In a previous paper [lj|, we calculated the charged-Higgs con- 
tributions to the mass splitting AM Bd and the decay rate E{B — > X s ^) at the NLO 
level, and found the strong constraints on free parameters X tt , A&& and Mjj from 
the well measured AM^ and E{B — > X s j). Two typical choices of (\ a , \bb) and 
the corresponding constraint on Mh obtained from the measured branching ratio 
of B — > X s 7 decay are 

III- A: (X tt ,Xbb) = (0.5,1), M^>150GeV (46) 

as shown in Fig. 9 of Ref.|13j|; and 

III - B : (X tt , X bb ) = (0.5, 22), 226 < M H < 285GeV. (47) 

For the first case, the new physics contribution to B — > X s ~f is very small and 
become negligible for M H > 250 GeV. For the second case (in [13j . it was denoted 
as the case C), the new physics contribution can be rather large, the sign of the 
dominant Wilson coefficient [m^] changed its sign from negative to positive 
due to the inclusion of the charged-Higgs penguin contributions. In this paper, we 
denote these two typical the model III-A and III-B, respectively. 



D. NLO Wilson coefficients in the general 2HDM's 

The new physics contributions to the quark level b — > s/d'y transition from the charged 
Higgs penguins manifest themselves from the correction to the Wilson coefficients at the 
matching scale Mw- In Ref . |37l| . the authors calculated the NLO QCD corrections to 
the B — > Xs'-y decay in the conventional models I and II. In Ref. l|> we extended their 



work to the case of model III. Here we firstly present the Wilson coefficients at the energy 
scales M\y and /i = 0(rrib) in a general 2HDM and then calculate the branching ratios, 
CP and isospin asymmetries, and the U-spin symmetry breaking of the exclusive decays 
B — > K*^ and B — > /ry in the following sections. 

Note that the CMM basis was used in Refs.[3, |3?|, the Wilson coefficients CfJ^ 
there are indeed the Wilson coefficients Zj(/i) in this paper. For the exclusive decays 
B — > V"f and to the first order in a s , only the NLO expression for C-j{fj) has to be used 
while the leading order values are sufficient for the other Wilson coefficients appeared in 
a?<y 7 ) in Eq.©. Therefore, only C 7 (/i), CfV) and zf\^) in Eq.© are affected 
by the charged-Higgs penguin contributions, while all other Wilson coefficients for i = 
1, . . . , 6 remain the same as in the SM. Since Z° 8 (/i) = C^/i), and C^(fi) = Z^(fi) j26{, 
so we here use the terms Cj S and C] for convenience. 



5 In this paper, the term model III always means the scenario of a general model III as presented in 
Ref.|4l|. In such model III 41], only the couplings Xtt and A^t, remain non-zero, and only the charged 
Higgs boson penguin dia gram pr ovide a new physics contribution to b — > 57 decay at one loop level. 
For more details see Refs.ll3, 411. 
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The new physics part of the LO Wilson coefficients Cj S at the matching energy scale 
Myy take the form, 

C° (M w ) = - \\Y\ 2 A{y t ) + (XY*)B(y t ) , (48) 



6' 

Cl NP {M w ) = -±\Y\ 2 D(y t ) + (XY*)E(y t ) , (49) 



where yt = mf/Mjj, and the functions A(x) and D(x) have been given in Eqs. fJlOJllj) . 
while 

3y - 5y 2 2y - 3y 2 

B{y) = : g[yi (50) 

The new physics parts of the NLO Wilson coefficients C^ s at the matching scale /iw 
can be written as 

Cl NP (M w ) = \Y\ 2 Cl YY {M w ) + (XY*) Cl XY (M w ), (52) 
Cl NP (M w ) = \Y\ 2 Cl YY (M w ) + (XY*) Cl XY (M w ) , (53) 



with 



Cl YY (M w ) = W i)YY + M i)YY In [y t ] + T i>YY Mn [x t ] - | 
Cl XY {M w ) = Wi, XY + M itXY ln[y t } + T i>XY (ln[x t ) - ^ J 



(54) 
(55) 



The explicit expressions of functions Wjj, Mjj and Tjj (i = 7, 8 and j = FY, J7 ) can 
be found in Refs. |37( or in Appendix B. The TV, terms appear when expres sing m t [Mw) 
in terms of the pole mass m t in the corresponding lowest order coefficients [371 ]. 

At low energy scale /i = 0(mb), the Wilson coefficients (7°' (//) and Cg(/i) after the 
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inclusion of new physics contributions can be written as 
C 7 V) = i? S [C 7jSM (Mw)+C° 7tNP (M w 



+ jj (^i - r,!) [C 8 ° SM (M W ) + C 8 ° iJVP (M w )] + ^ V ai , (56) 

i=i 

8 

CSO*) = [C° 8 , SM (M W ) + Cl NP (M w )} + ft ^ , (57) 

i=l 

C 7 V&) = ryi [C 7 1 iSM (M iy ) + C 7 1 jiVP (M w )] 



+ 3 (tj* -17* J [CS, SM (M W ) +C 8 1 )JVP (M W )] 

/ 297664 16 7164416 u 256868 37 6698884 39 
+ I 7723 — 7723 -| 7723 7723 

V 14283 ' 357075 ' 14283 ' 357075 ' 
■ [Cl SM (M w ) + Cl NP (M w )] 

+ W - V§ ) \- C ^m(M w ) + Cl NP {M w )] 

8 

+ J2 hvE(x t ) + fi + giT])^} V ai , (58) 

i=l 

where the "magic numbers" are listed in Table |H] 

In the conventional model I and II, the general Yukawa couplings X and Y are real 
and given by 

X = -cot p, Y = cot(3 (Model I), (59) 
X = tan/?, Y = cot{3 (Model II) . (60) 

In the Model III where only the couplings Xu and A&& are non-zero, the relation between 
the couplings (X,Y) and (Ait,A&&) is also simple 

X = -X bb , Y = \ tt (Model III). (61) 

Now we are ready to calculate the numerical results for the B — > decay in the 
general 2HDM's with the inclusion of NLO QCD corrections. 



III. B -► iT* 7 DECAY 

For the numerical calculations, unless otherwise specified, we use the central values 
of the input parameters as listed in Table and consider the uncertainties of those 
parameters as given explicitly in Table IIHI 

From Eqs.()30|) and (|38|h the decay amplitude and branching ratio for B — > K*^ decay 
can be written as 

A(B - K*j) = ^R K * (K*j\Q r \B) , (62) 



V2 



G 2 F am\m\ f _ m|* \ 3 , T) , , , , , , ., 
32tt 4 1 m\ 



B(B -> K*j) = T B ~*~ n J b ( 1 - ^ ) \Rk*\ Wk\\ (63) 
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with 



Rk* = V* s V ub [a u 7 (K*~f) + a: nn (K*j)] + V c * s V cb [a£(iT 7 ) + < nn (K* 7 )] ■ (64) 
The CP asymmetry of B —* K*-y can also be defined as [Til , ^( 

v ; r(_B — > x*7) + t(b — > if 7) v ; 

Another physical observable for B — > V7 decay is the isospin symmetry breaking in 
the K*^ — K* or p ± — p° system. Since the branching ratios of both B — > K*^^ and 
B -> if 7 decays have been measured, the study of the isospin breaking in B — > ^7 



decays becomes very interesting now |27l. l28|. Following Ref . |27l| . the breaking of isospin 

*o 

symmetry in the K*~ — K system can be defined as 

s , rB (^r%)- B (^A- 7 ) ^ 
t? t S(5 -> if 7) + £(£ -> if*- 7 ) 

where r/ T = t b +/t b o, and the CP-averaged branching ratios are understood. 

By using the world averages as given in Eq.(j2} and the ratio r B + /t b o = 1.083 ± 0.017 
(3jJ, we find numerically that 

A -(iT 7 r p = (3.9 ± 4.8)%, (67) 

where the errors from the two measured branching ratios and the ratio r B + / t b o have been 
added in quadrature. The measured value of isospin symmetry breaking is indeed small 
as expected previously 0, Any new physics contribution producing large isospin 
breaking for B — > if* 7 decays will be strongly constrained by this measurement. 



A. Branching ratios and CP asymmetries 

By using the formulas as given in Eqs. ()32|) and (}3*3*j) and the central values of input 
parameters in Table ITTT1 we find the SM predictions for a 7 {K*^) and a ann (K*j) at the 
low energy scale \i = m b , 

a%(K*j) = -0.3212 + 00113 -0.1407 - 0.0683z 0.0330 - 0.0002i 

C^ SM (mf,) ^^7,SAf T 1 —contribution T 11 —contribution 

= -0.4177 - 0.0685i, (68) 
a c 7 (K*j) = -0.3212+^0113 -0.0802 - 0.013H r 0.0161 - 0.0120i 

SM (mt,) ^^7,SM T 1 —contribution T 11 —contribution 

= -0.4063 - 0.0251i, (69) 

a u ann {K 7) = < nn (K 7) = -0.0092, (70) 
a u ann (K*-j) = 0.1933, (71) 
a c ann (K*-j) = 0.0046. (72) 
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It is easy to see that (a) the type-I contribution is about 4 times larger than the type- 
II contribution, and (b) only the weak annihilation factor d£ nn (K*~ 7) contributes to 
the decay B — > K*^ effectively, since for b — > u transition the power suppression is 
compensated by the large Wilson coefficient C\ and the occurrence of annihilation at tree 
level. 

The corresponding NLO SM predictions for branching ratio BiB — > K*j) are 

B(B^T°j) SM = [3.35^;^(F^)^(^)^o 7 (Ab)±0.20K)] xlO- 5 

= (3.35±i;H) x lO" 5 , (73) 

B(S^^*- 7 ) SM = [3.25^S(Fx.)18^(A*)lg:S(AB) ± 0.20K)] x 10~ 5 

= (3.25iiS) x lO" 5 , (74) 

where the four major errors have been added in quadrature. The uncertainty of the form 
factor Fx* dominate the theoretical error, and the remaining errors from other input 
parameters are negligibly small. Although the central values of the SM predictions for 
the decay rates are smaller than the world average as given in Eq.(J2J), but they are in 
good agreement within la theoretical error. The effect of annihilation contribution on 
the decay rates is less than 5% numerically. 

If we use F K * = 0.38 ±0.06 obtained from the light-cone sum rule (LCSR) |32| instead 
of Fk* = 0.25 ± 0.06 in numerical calculation, we find 

£(£> k*°j) sm = (7.27±£jj?) x lO" 5 , (75) 
B(B -> K*-j) 3M = (7.3li|g) x 10" 5 . (76) 

Here the central values are much larger than the measured values as given in Eq.(j2J), 
but still agree with the data within 2a errors because of the large theoretical error. For 
the purpose to study the new physics contributions to the exclusive decays B — ► V7, 
one prefers a better agreement between the SM predictions and the high precision data. 
Therefore, we will use Fx* = 0.25 ± 0.06 in this paper, unless otherwise specified. 
For the model I, the theoretical predictions for branching ratios are 

B(fl-irV = (3.35±i:H) x lO" 5 , (77) 
B(B -> K*-^) 1 = (3.25i{;^) x 10" 5 , (78) 

for tan/? = 4 and Mh = 200 GeV. The Mh dependence of the branching ratio B(B — > 
K*^f) is very weak: it will change by less than 2% in the range of 200 < Mn < 600 GeV. 

Figffl shows the tan j3 dependence of the branching ratio BiB — > K*^^) in model I 
for Mh = 200 GeV. The dots and solid line shows the central value of the NLO SM 
and model I prediction, respectively. The region between two dashed lines shows the 
NLO SM prediction with error as given in Eq. (|74p . The shaded band shows the data: 
B(B -> K*--f) exp = (4.18 ± 0.32) x 10~ 5 . From this figure, one can see that (a) the NLO 
SM prediction agree with the data within la error; and (b) the new physics contribution 
in model I is negligibly small for tan/5 > 1, while a tan/3 < 0.5 is also strongly disfavored. 
For B —>■ K* 0/ -f decay mode, we have the same conclusion. 
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In the popular model II, the numerical results for a^{K*^f) at the low energy scale 
/i = nib are, 



a^(K*-i) 11 = -0.3100 - ^06523, - 0.1436 - 0.0724z +0.0481 - 0.0003i 

CV,SM( m f>) ^C7,jvp T 1 —contribution T 1 1 —contribution 

= -0.4707 - 0.0728i, (79) 
a c 7 (K*-i) 11 = -0.3100-^0652 -0.0831 - 0.0172j -0.0265 - 0.0182i 



C7,SM{rrib) ^C-7,NP T 1 —contribution T 1 1 —contribution 

= -0.4848 - 0.0354i, (80) 

for tan {3 = 4 and Mh = 300 GeV. The second terms in above two equations are the new 
physics corrections to the NLO Wilson coefficient CV^M^b); the hard vertex and hard 
spectator contributions are also changed slightly because of the variations of Z®(fi) and 
C®([ih) after including the charged-Higgs contributions. The total new physics contribu- 
tion to ay in model II is around 10% for tan (3 = 4 and Mh = 300 GeV. The annihilation 
parts remain unchanged. 

For the model II, the theoretical predictions for branching ratios are 

B{B^T°>y)» = [4.54^(F^)^(^)^(Ab)±0.22K)] xlO- 5 

= (4.54^;1) x lO" 5 , (81) 

B{B-+K-i)» = [4.47l?;H(i^.)18:P)(^)^(Afl)±0.23(m c )] x KT 5 

= (4.47+^6) x io-5 ; (82) 

for tan (3 = 4 and M H = 300 GeV. 

Fig. El shows the Mh dependence of the branching ratio B — > K*^ in model II for 
tan (3 = 4. The dot-dashed and solid curve shows the central value of the NLO model II 
prediction for the branching ratio B{B — > K 7) and B[B — > K*~j), respectively. Other 
band or lines show the same thing as in Fig. ^ 

Fig. El shows the tan (3 dependence of the branching ratio B — > K*^ in model II for 
Mh = 300 GeV. The dot-dashed and solid curve shows the central value of the NLO model 
II prediction for the branching ratio B(B — > K 7) and B[B — > K^'j), respectively. 
Other band or lines show the same thing as in Fig. ^ 

It is easy to see from FigO that a charged Higgs boson with a mass around 200 
GeV is still allowed by the measured branching ratio of the exclusive B —>■ K* / y decay, 
which is weaker than the lower bound of Mh ^ 300 GeV obtained from the data of 
the inclusive B — > Xg'j decay. This is consistent with general expectation. The key 
point here is the large uncertainty of the non-perturbative form factor Fx*. If we use 
Fr* — 0.38 ± 0.06 and keep all other input parameters unchanged, we get a much strong 
lower limit on Mh, as can be seen from Fig 01 where the solid and dot-dashed curves show 
the NLO model II prediction for B{B — > K 7) and B{B — > K^'j), respectively. The 
dots line and the band between two dashed lines show the corresponding SM prediction 
of B{B -> K*~i) = (7.3lj|;p) for F K * = 0.38 ± 0.06. 

Now we study the model III. According to previous studies in Ref . jljj] . we got to know 
that the charged Higgs penguins can provide a significant contribution to the dominant 
Wilson coefficient CV(/x) and changed its sign from negative to positive. Of course, the 
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size of the new physics contributions is strongly constraine d by the measured branching 
ratio of the inclusive B — > Xgj, as investigated in detail in [13j | . 

For the model III-A, i.e. (X tt , \ bb ) = (0.5, 1), the new physics contributions are small, 
the numerical results for a,j(K*^f) at the low energy scale \l = nib are, 



a^(K*-f) m ~ A = -0.3100 +^029g - 0.1394 - 0.0664i +0.0336 - 0_0002i 

C7,sm(ti6) AC7,jvp T 1 —contribution T 1 1 —contribution 

= -0.3859 - 0.0666i, (83) 
a c 7 (K*7) m ~ A = -0.3100+^0293 j-0.0788- 0_0112i -0.0155 - 0.01202 

CV,SM("ii>) ^C7,JVP T 1 —contribution T 11 —contribution 

= -0.3744 - 0.023H, (84) 

for Mh = 300 GeV. The total new physics contribution to a 7 in model III-A is also around 

10% in magnitude for Mh = 300 GeV, but in the opposite direction of that in model II. 

The annihilation parts also remain unchanged. 

Fig. 03 shows the Mh dependence of the branching ratio B — > K*^y in model III-A. 

The dot-dashed and solid curve shows the central value of the NLO model III prediction 

*o 

for the branching ratio B[B — > K 7) and B{B — > K*^^), respectively. Other band or 
lines show the same thing as in Fig. ^ The new physics contribution here is small are 
consistent with the SM prediction within la error. Numerically, we have 

B (B -> Z* 7 ) m - A = (2.871J-J?) x 10" 5 , (85) 
B(B -> K*- 7 ) m - A = (2.75iJ;^) x 10" 5 , (86) 

for Mh = 300 GeV, where the four major errors as in Eqs. (j8H) and (J8*2j) have been added 
in quadrature. 

For the model III-B, i.e. (At*, Aw,) = (0.5,22), the new physics contributions are large, 
the numerical results for a?(K*j) at the low energy scale fi = nib are, 

a^(K*7) m " B = -0.3100+^8485 r 0.1049 - 0.0174i +0.0752 - 0.0003i 

Ci,SM{ m b) AC-7,NP T 1 —contribution T 11 —contribution 

= 0.5088 -0.01772, (87) 
a c 7 (K*-f) m - B = -0.3100+^8485 r 0.0443 + 0.0379? +0.0005 - 0.0182i 

Cr,SM(mb) ^-C-j t NP T 1 —contribution T 11 —contribution 

= 0.4947 + 0.0196i, (88) 

for Mh = 250 GeV. The second terms in above two equations are the new physics cor- 
rections to the NLO Wilson coefficient C7 iSM (m&), which is large and positive and makes 
the a 7 (K* , ~f) to be positive also. The hard vertex and hard spectator contributions are 
also changed moderately, but has only small effects on the branching ratios. 
For model III-B, the theoretical predictions for branching ratios are 

£(B - 7TV n ~ B = [4.23™(F^)^(^)^(A B )±0.21K)] x KT 5 

= (4.23l£g) x 10- 5 , (89) 

B(B - K^r-* = [5.07i^(F^)^)(^)i^(AB)±0.24(m c )] x KT 5 

= (5.07™) x 10- 5 , (90) 
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for M H = 250 GeV. 

Fig. El shows the M# dependence of the branching ratio B — ► if*7 in model III-B. The 
dot-dashed and solid curve shows the central value of the NLO model III-B prediction for 
the branching ratio B[B — > K 7) and B{B — > K*~^f), respectively. Other band or lines 
show the same thing as in Fig. [SJ 

If we add the theoretical errors as given in Eqs. 1)89190)) with the corresponding experi- 
mental errors in Eqs.(j2J) in quadrature and treat them as the total la error, we then read 
off the allowed regions of M# from Fig El 

218 <M H < 293GeV, and M H > 1670GeV, (91) 

allowed by the measured B{B — > K 7), and 

226 <M H < 315GeV, and M H > 1490GeV, (92) 

allowed by the measured B(B — > K*~^). These constarints on Mh are well consistent 
with those obtained from the inclusive B — > X s 7 decays as given in Ref.j3|. Of course, 
the large theoretical error is dominated by the uncertainty of the form factor Fk* here. 

For the exclusive B —>■ K* / y decay, the theoretical prediction for the CP symmetry 
Acp as defined in Eq.(|65|) is very small: 

\A CP (B ^K* 7 )l < 1% (93) 

in the SM and all three types of the 2HDM's considered here, which is consistent with 
the measurements as reported by BaBar ^| and Belle Collaboration [l6| : 

A C p(B -> iT 7 ) = [-0.17, +0.082], (94) 
A C p(B -> K*-f) = -0.001 ±0.044 ±0.008. (95) 



B. Isospin Symmetry 

As can be seen in last subsection, the large uncertainty of the form factor Fk* domi- 
nates the total error of the theoretical prediction of the branching ratios. For the isospin 
symmetry breaking of B — > K*^ system, however, its dependence on the form factor Fk* 
largely cancelled in the ratio. From Eqs. (|38l31|) . the isospin symmetry breaking A __(K*~f) 
as defined in Eq . ()66j) can also be written as 



A -(^*7) = 



Vr\ 




2 _ 


R K *- 


12 


Vr 


1%*° 


l 2 ± 


\Rk*° 1 


2 



(96) 



where Rk* have been given in Eq.(|64|). In our approximation, the isospin breaking is 
generated by weak annihilation contributions, and has a residue sensitivity to the form 
factors Fy induced through the F v dependence of b v and d v functions as defined in 
Eqs. (|36l37|) . Since = V* s V u b is about two orders smaller than = V* s V c b, the function 
Rk* is largely determined by a^{K*^). 
In the SM, we have numerically 

A -(iT7) SM = [5.6il;[(F^)^:i(^)ii1(AB)±0.1(m c )] x lO" 2 

= (5.61*5) x 10- 2 , (97) 
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where the errors added in quadrature. The dominant error comes from the uncertainty of 
the low energy scale l/2m& < fi < 2m;,. The SM prediction agrees well with the measured 
value of A e */(K*-f) = (3.9 ± 4.8)%. 

In the general two-Higgs-doublet models, by assuming tan (3 = 4 and Mh = 250 GeV, 
we find numerically 



where the errors induced by the uncertainties of u.,Fk* } \b and m c have been added in 
quadrature, and the uncertainty of fi dominates the total theoretical error. 

Fig. [7| shows the [i dependence of the isospin symmetry breaking A _ {K*~f) in the 
general 2HDM's for tan f3 = 4 and Mh = 250 GeV. Two coincided dot-dashed curves show 
the SM and model I prediction , the dash and dots curve show the model II and model 
III-A prediction respectively, and the solid curve refers to the model III-B prediction. The 
error bar shows the data A _(K*j) exp = (3.9 ± 4.8)%. 

Fig. IS] shows the Mh dependence of the isospin symmetry breaking A _ (K*j) in the 
general 2HDM's for tan (3 = 4 and \i = GeV. Two coincided dot-dashed curves show 
the SM and model I predictions. The dashed, dots solid curve show the model II, III-A 
and III-B prediction, respectively. The error bar shows the data as in Figd 

From above two figures, one can see that only the theoretical prediction of the model 
III-B is rather different from that of the SM and looks like deviating from the data. But 
the regions M H < 200 GeV and 300 < M H < 1500 GeV have been excluded by the data of 
inclusive B — > [l3| and by the constraint as illustrated in FigEl The main reason for 
the great changes of the solid curve in Fig|H]is the strong cancellation between the negative 
CV,sA/( m &) and its positive new physics counterpart as illustrated clearly in FigEl where 
the solid curve shows the summation of the SM and new physics contributions to the 
dominant Wilson coefficient CV, i.e., CV(mf,) = CV )SM (m&) + ACV jJVP (mb). When CV(m&) 
approaches zero, the summation of other "originally small" parts (such as the T/ , T/ 7 
and a p ann contributions ) becomes important and leads to an abnormally large isospin 
breaking. The short-dashed and dot-dashed curves in FigEl shows the absolute value 
of R-j^*o and Rk*-, respectively. The isospin breaking is proportional to the difference 
of their squares. At the close region of the crossing point of R-^*o and Rk*-, the ratio 
A _(i^*7) can be large and changes the sign. But as mentioned previously, this region 
around Mh = 500 GeV has been excluded by the data of branching ratios from both the 
inclusive and exclusive radiative B meson decays. 

In the region of Mh ~ 250GeV, the model III-B is disfavered by the measured value 
of isospin breaking Ao_ (K*^) as can be seen from Figs. [7| and |H1 But taking the sizable 
experimental and theoretical uncertainties into account, the theoretical prediction of the 
model III-B is still compatible with the data within 2a errors. In another word, the 
positive CV(m& ) is disfavored but can not be excluded by the present data. 



Ao-(iT 7 ) 
Ao-^7) 1 
A -(^7) m -' 
Ao-tiTT) 11 " 



(+5.7j£?) x 10- 2 
(+4.6+1;!) x io~ 2 
(+6.2^;o 7 ) x 10- 2 

(-5.1135) x io- 2 



(98) 
(99) 
(100) 
(101) 
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IV. B -» py 



When compared with B — > K*^ decay, the B — > p7 decay mode is particularly in- 
teresting in search for new physics beyond the SM, because of the suppression of b — > d 
transitions in the SM and the simultaneous chirality suppression. For B — > fry decay, 
we generally know that: (a) both a" and a 7 contribute effectively since A^ and Af are 
comparable in magnitude; (b) the branching ratios of B — > p7 are suppressed with re- 
spect to -B — > fC*7 by roughly a factor of | V^^/ Vts 1 2 ~ 4 x 10~ 2 ; (c) the CP asymmetry 
Acp{B — > fry) is generally at 10% level, and may be observed in B factory experiments; 
(d) the new physics may provide significant contribution to the observables of B — > p7 
decay; and (e) only the experimental upper limits on the branching ratios of B —>■ p7 are 
available now. 



A. Branching ratios and CP asymmetries 

From Eq. (j38)l . the branching ratios of B —>■ pj decays can be written as 



G \am\m\ ( _ rrf x " 
32tt 4 1 mV 



B(B - fry) = t b ~*ZT > ( 1 " ~T ) I*pI ^p| 2 , (102) 



with 



^p = [«r(P7) + Cn(P7)] + V; d V cb [a c 7 (p 7 ) + < nn (p 7 )] . (103) 

By using the formulas as given in Eqs. (|3*2*|) and (f3"4*|) and the central values of the input 
parameters in Table ITTT1 we find the SM predictions for a 7 (p7) and a ann (fry) at the low 
energy scale p = m b , 

a%{pj) = -0.3212 + 00113 r 0.1407 - 0.0683i 0.0343 - 0.0003i 

^7SAf( m f>) ^^7,SM T 1 —contribution T 11 —contribution 

= -0.4164 -0.0686z, (104) 
a c 7 {fry) = -0.3212 +^0113, - 0.0802 - 0.013b; r 0.0166 - 0.0143i 

SM (rn b ) ACj SM T 1 —contribution T 11 —contribution 

= -0.4067-0.02742, (105) 
Cn(P°7) = -0.0032, a c ann (p°j) = -0.0127, 

<nn(P~l) = 0-1883, < nn (p- 7 ) = 0.0032. (106) 

Here the values of weak annihilation factors are slightly different from those for B — > K*^ 
decay, and only the T 11 contributions to a 7 (p7) are different from those to dj(K*^) because 
of the small differences of the Hi functions between two decay modes as can be seen in 
Appendix A. 

The corresponding NLO SM predictions for branching ratio B(B — > fry) are 

B{B -> p° 7 ) SM = [0.91 ± 0.29(j)t° 2 i(F p ) ± 0.17(p)i^(A B ) ± 0.10(m c )] x lO" 6 

= (0.911°;g) x 10- 6 , (107) 

B(B - p" 7 ) 5M = [2.03 ± 0.34( 7 )^:f 7 (F p ) ± 0.31(p)^(A B ) ± 0.12(m c )] x 10~ 6 

= (2.03±g;||) x 10~ 6 , (108) 
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where the individual errors have been added in quadrature. The uncertainties of the CKM 
angle 7 ( here we take 7 = (60 ± 20)° in the calculation) and the form factor F p dominate 
the total error, and the remaining errors from other input parameters are negligibly small. 

The central values and theoretical uncertainties of the branching ratios B(B — > /ry) in 
the SM and the general 2HDM's are all listed in Table IIVI The SM prediction is well 
consistent with the experimental upper limits within la error. The predictions of the 
model I , II and III-A are also compatible with the data within errors, as illustrated in 
FigEl and FigE] for B — > p°7 and B — > p~7, respectively. For model III-B, however, 
the branching ratios can be changed significantly when the charged Higgs boson is light 
or heavy, as illustrated by the solid curves in FigES and FiglTTl From the experimental 
upper bounds on B(B — > p^y) as given in Table HI we find the lower limit on Mh 

M h > 206GeV (109) 

when the 2a theoretical errors are also taken into account. This lower bound is compatible 
with those obtained from the measured BiB — > K*j) as given in Eqs. (j91|) and (|92|) and 
from the inclusive B — > X s ^ decay 

The CP asymmetry of B — > /ry decays is defined in the same way as for B — > K*^ 
decays in Eq. (j65J) . Using the input parameters as listed in Table UTTl one finds the NLO 
SM predictions, 

[SAtlM ± 1.9(R b ) ± 0.8(A B )1?;?K) ± 0.4(F p )+?i( 7 )] x lO" 2 
(8.41^) x 10- 2 , (110) 
[1O.4JK00 ± 2A(R b )t 3 (X B ) ± 0.8(m c ) ± 0.1(F p )l?j( 7 )] x 10" 2 
(10.4±S5) x 10- 2 , (111) 

where the errors have been added in quadrature. The CP asymmetry of B —>■ pj is large 
is size and depends sensitively on the variations of the scale p and R b = \/p 2 + fj 2 . If we 
consider the whole range of 0° < 7 < 180° instead of 7 = (60 ± 20)° preferred by the 
global fit result [31] , the CP asymmetry Acp(B — > pj) also shows a strong dependence 
on the angle 7 as illustrated by FigH2]for B ± — > p ± 7 (solid curve) and B — > p°j (dashed 
curve) decays. 

The numerical values of CP asymmetries in the SM and the general 2HDM's are also 
listed in Table IIVI The theoretical predictions of the SM and the model I, II and III-A 
are all compatible, around +10%. The CP asymmetry Acp{B — > pry) in model III-B, 
however, is comparable in size with the SM prediction, but has an opposite sign, as 
shown in Fig |13| where the upper and lower three curves show the theoretical predictions 
for p = mb/ 2 (dashed curves), m h (solid curves) and 2rrib (dots curves) in the SM and the 
model III-B, respectively. For the B — > p°7 decay mode, we have the similar conclusion. 
This feature may be served as a good observable to distinguish the model III-B (or a 
positive C 7 (rrib)) with the SM (a negative C-jirrib))- 



Acp(p°i) sm = 



B. Isospin and U-spin symmetries 

According to currently available data, the SU(2) isospin symmetry of the strong inter- 
action is a very good symmetry with a breaking no more than 5%. The U-spin symmetry, 
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TABLE IV: The NLO theoretical predictions for branching ratios and CP asymmetries in the 
SM and the models I, II, III-A and III-B, assuming tan (3 = 4 and Mg = 250GeV. The errors 
induced by the uncertainties of six input parameters (p, Rj,, \B,m c , F p and 7) are taken into 
account. Individual errors are added in quadrature. 



Decays 


SM 


model I 


model II 


model III-A 


model III-B 


0(5° _> p° 7 )(10- 6 ) 


u ' yi -0.40 


n qn+ - 41 

u - JU -0.39 


1 sn+ - 59 


u -' D -0.33 


1 n7 +o.50 


B(B-^p- 7 )(10- 6 ) 


2 +0 - 8 


2 +0 ' 8 


2 q +L2 


1.718:5 


2 4 +L4 


Acp{^^ A)(%) 




8.5+f-| 


7-0±l;g 




7 q+3.6 
~ ' • / -4.6 


Acp(B- - p- 7 )(%) 


10.4l|° 


10.511;° 


8.71^ 




tx+4.2 

-»-t>-5.o 



the SU(3) flavor symmetry of the strong interaction under exchanges of the down and 
strange quarks, however, may have a breaking around 20% (i.e., ~ (F^/F^ — 1) ) as 
frequently used in the study of B — > Kn decays. For the exclusive B — ► K*j decays, 
the isospin breaking derived from the measured branching ratios is indeed around 5% as 
given in Eq. (|67l) . For B — > p7 decays, no measurements are available now. 

As in Refs. |24l. l2r| , we also define the isospin symmetry breaking of B — > p7 decays as 
the form of 



A(P7) = 2 



r(B+ -> p+ 7 ) r( J e+ -> p+ 7 ) 



+ 



2r(5 -> p° 7 ) 2r(5° -> p° 7 ) 



(112) 



Using the central values of input parameters as listed in Table ITTT1 and assuming tanfi 
Mh = 250 GeV, we find numerically that 



A(p 7 ) 



(0.91^) 


x 10" 


-2 


in 


SM, 






(0.9 'I-) 


x 10~ 


-2 


in 


model 


I- 




(0.4™) 


x 10" 


-2 


in 


model 


II, 




(i-3±!!;l) 


x 10" 


-2 


in 


model 


III 


-A, 


(4.91-°) 


x 10" 


-2 


in 


model 


III 


-B, 



(113) 



where the errors from uncertainties of input parameters have been added in quadrature. 
The largest theoretical uncertainty comes from the the CKM angle 7. 

In Figim we show the angle 7 dependence of the isospin breaking A(p 7 ) in the SM 
and the considered 2HDM's for tan/5 = 4, M H = 250 GeV and 0° < 7 < 180°. It is easy 
to see from FigE] that (a) except for the model III-B, the isospin breaking in the SM 
and other 2HDM's have the similar 7 dependence; (b) all theoretical predictions become 
almost identical and very small in magnitude for 7 ~ 55° (the value preferred by the global 
fit results), and the smallness of A(pj) is also consistent with the general expectation and 
other measurements; (c) the theoretical predictions in the SM and model III-B have a 
very different 7 dependence, and have the opposite sign for small or large values of the 
CKM angle 7. 

The U-spin symmetry is another interestin g ob servable for B — > (iT*,p) 7 decays, and 
has been studied for example in Refs. [2J, |26|, |45j . In the limit of U-spin symmetry, the 
quantity 



AU(K*,p) = AB{B -> iT 7 ) + AB{B -> p 7 ) = 



(114) 
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with 



AB(B -> K*j) = B(B H 
AB{B -> p 7 ) = fi(5 H 



p+ 7 ) - B(B- 



(115) 
(116) 



should be satisfied. Using the central values of input parameters, we find the SM predic- 
tion of AU(K*, p) 



AB(B -> K*i) 
AB(B -> p 7 ) 



-3.7 x 10~ 7 , 
+4.4 x 10~ 7 , 



(117) 
(118) 



where we have chosen 7 = 90° which maximizes the effects. The two parts have opposite 
sign and cancels to a large extent, leaving a small U-spin breaking 



AU(K*,p) = 0.7 x 10" 



(119) 



in the SM, which is only about 8% of BiB — > p°j). In the general 2HDM's, we find the 
numerical results 



AU(K*,p) 



0.7 x 10~ 7 in model I, 

0.9 x 10~ 7 in model II, 

0.6 x 10~ 7 in model III - A, 

-1.5 x 10~ 7 in model III - B, 



(120) 



for tan/3 = 4, Mh = 250 GeV and 7 = 90°. The new physics contributions in the 
conventional model I , II and the model III-A have little effect on the size of U-spin 
symmetry breaking. In model III-B, although AU(K*,p) becomes negative, but it is still 
small in magnitude. 



V. CONCLUSIONS 

By employing the QCD factorization approach for the exclusive B — > V7 decays 
as proposed in Refs. [22I I23L 24]|, we calculated the NLO new physics contributions to 
the branching ratios, CP asymmetries, isospin symmetry breaking and U-spin symme- 
try breaking of the exclusive radiative decays B — > K*j and B —>■ fry, induced by the 
charged Higgs penguin diagrams appeared in the general two-Higgs-doublet models in- 
cluding the conventional model I and II, as well as two-typical cases of model III. The 
NLO new physics contributions are included through their corrections to the NLO Wilson 
coefficients C 7 (M W ) and C S (M W ) at the matching scale M w . 

In section |Hj we gave a brief review about the effective Hamiltonian and the calcu- 
lation of the exclusive B — > U7 (V = K*,p) decays at next-to-leading order in QCD 
factorization, presented the relevant formulas for the calculation of Wilson coefficients 
and physical observables in the SM and the general two-Higgs-doublet models. 

In section ITTT1 and ITVl we calculated the NLO new physics contributions to the branch- 
ing ratios and other observables of B — > K*^ and B — ► p7 decays in the general 2HDM's, 
compared the theoretical predictions with those currently available experimental mea- 
surements, and found the following points: 
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The new physics corrections to the physical observables under consideration in this 
paper are generally small in the model I and model III-A, moderate in model II, but 
large in model III-B. And therefore the theoretical predictions in the SM, model I 
and III-A are always in good agreement with the corresponding data. 

For model II, a lower bound on the mass M# can be obtained from the measured 
branching ratios of B — > K*^ decays: 

M H > 200 or 300GeV, (121) 

if one uses Fk* = 0.25 ± 0.06 or Fk* = 0.38 ± 0.06 in the calculation, as illustrated 
by FigJ2]and FigEJ From Fig|31 a lower limit of tan/5 > 0.5 can also be obtained 
from the data. 

In the model III-B, the new physics contributions to C 7 $(My/) are larger than its 
SM counterparts in size and change the sign of the dominant Wilson coefficient 
C 7 (m b ) from negative to positive, as given in Eqs. (j8Tj) and (J581) . 

In the model III-B, the ranges of 

226 <M H < 293GeV, and M H > 1670GeV, (122) 

are still allowed by the measured B(B — > K*^y) as given in Eq.Q. From the 
experimental upper bounds on B{B — > fry), we find the lower limit on Mh 

M H > 206GeV (123) 

when the 2a theoretical errors are also taken into account. Above limits on Mh are 
comparable with those obtained from the inclusive B — > X s ^ decay 

The model III-B prediction for the isospin symmetry breaking of B —>■ K*^ decay 
is A _(K*'j) = (— 5.6l|g)%, which is small in size but has an opposite sign with 
the measured value, as illustrated in Figs. [7| and |H1 A positive CV(m{,) is therefore 
disfavered by the measured value of Aq^(K*j) = (3.9 ±4.8)%, but still can not be 
excluded if we take the large theoretical and experimental errors into account. 

The theoretical predictions for CP asymmetry Acp{B — > K*^) is always less than 
one percent in magnitude in the SM and all three types of the general 2HDM's 
considered here. For B —*■ fry decay, however, its CP asymmetry can be as large 
as about 10% in size in the SM and all three types of 2HDM's and have a strong 
dependence on the variations of the scale /i = 0(rrib) and the CKM angle 7, as shown 
in Figs. and EH It is interesting to see from Fig^] that the CP asymmetry in 
model III-B has an opposite sign with the one in the SM. This feature may be used 
as a good observable to distinguish the model III-B (or a positive C 7 (m b )) with the 
SM (a negative C 7 (rrib))- 

For B — > fry decay, the isospin symmetry breaking is less than 10% in the region of 
7 = [40 ~ 70]° as preferred by the global fit result [3l[, but can be as large as 20 
to 40% in the regions of 7 < 10° and 7 > 120°, as can be seen clearly in Fig. 
The SM and model III-B predictions for isospin breaking have an opposite sign for 
small or large values of the CKM angle 7. 
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• The U-spin symmetry breaking AU(K*,p) in the SM and all 2HDM's considered 
here is generally small in size: ~ 10~ 7 . 
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APPENDIX A: d AND H( FUNCTIONS 

In this Appendix, the explicit expressions or numerical values of all Gi and Hj func- 
tions appeared in Eq. (jHU|) will be listed. For more details of these functions, one can see 
Ref. 12 g| and references therein. 



G 1 {z) 
G 2 (z) 
G 3 
G 4 (z c ) 

G R ( Zr ) 



G* 



52 lnA 

81 rrib 
104 



833 1. , , , , . . 10z7r 

\a(z) + b(z)] H . 

972 4 L w WJ 81 ' 



, ,, , a 833 3 r , . , . , , 20iir 
27 rrib 162 2 27 



27 rrib 

38 , u 
— In — + 

81 m b 

1568 , a 
In — + 

27 m b 



598 2tt 8 3 3,, lN 14ztt 

h — = + -X & a(l) + -6(1) H . 

81 V3 3 4 w 2 W 27 



761 
972 ~ 
14170 



81 

2855 



7r 4 1 5,. . 37i7r 

— = X b + -a(l) + -b(z c ) 

3x/3 9 8 V ; 4 v cj 81 

8vr 32 . , . 224«tt 

+ + y X b - 12a(l) + 246(1) + — , 



1156, u 
ln — H 

81 m b 486 



47T 

3^ 



9 



~a(l) + 116(1) + 9a(z c ) + I5b(z c ) - 

L (51 

8, a 11 2vr 2 2m 

-In— H 1 , 

3 3 9 3 



(Al) 
(A2) 
(A3) 
(A4) 
(A5) 

(A6) 

(A7) 
(A8) 



where 



a(l) 



dx dy dvxy\n[v + x(l — x)(l 
Jo Jo Jo 



Am 

4.0859 + —, 
9 



6(1) = ®° 
v ; 81 



4tt 632tt 2 
+ 



3^ ' 1215 15 

a( Zu ) = (-1.93 + 4.96?) x 10~ 5 , 

a(z c ) = 1.525 + 1.242i, 

b(z u ) = (1.11 + 0.28?) x 10" 5 , 

b( Zc ) = -0.0195 + 0.13182, 



d 2 lnT(x) 

dx 2 



v)(l — v + vy)] 



Am 



-0.1684, (A9) 

(A10) 

0.0316 + —, (All) 

(A12) 
(A13) 
(A14) 
(A15) 



where z q = m 2 /m 2 and the masses m q (q = u,c,b) as listed in Table ITTT1 have been used 
to obtain the numerical results. The explicit analytical expressions for a(z) and b(z) can 
be found for example in Ref.j26[. 
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For the HY functions, we have 



HY(z p 



H 



27T 2 f B f^X B 
9 FyVtlB 



dv h(v, z p )§y(v) 







Hj{z c 



v 



Hi 



--[HY(l)+H?(0)} 

Hl{z c ) - 1^(1), 
2^(1), 

-hX{z c ) + \hX0) = 

^ 2 f B f^X B ^ 



3 F v m B 

where the hard-scattering function h(u, z) is given by 



v i v i 
a{ + «2 + 



Az 

h(u, z) = — { Li 
u 



1 



u—iz+ie 



+ Li-2 



1 + 



u—iz+ie 



2 

u 



(A16) 

(A17) 
(A18) 

(A19) 
(A20) 
(A21) 

(A22) 



(A23) 



where Li 2 [x] is the dilogarithmic function, and the function h(u, z) is real for u < Az and 
develops an imaginary part for u > Az. The light-cone wave function takes the 

form of 



$±( v ) = 6v(l -v)\l + aY(/j)C 3 1 /2 (2v - 1) + a\ \/j)C 3 2 /2 \2v - 1) + 
where C? /2 (z) = 3x, C 2 3/2 (x) = |(5a; 2 - 1). 



(A24) 



APPENDIX B: NLO COEFFICIENTS AT /x = My^ IN GENERAL 2HDM'S 



For the completeness, we list here the expressions of the NLO functions Wij,Mij 

7, 8 and j = YY, XY) at the matching scale fiw = Mw in the general 
two-Higgs-doublet models. For more details see Ref. |37j]. 
The NLO functions proportional to the term \Y\ 2 are 



and Tij (i 
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8y 3 
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The NLO functions proportional to the term (XY*) are 
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FIG. 1: Plot of the branching ratio B(B -» K*-j) vs tan (3 in model I for M H = 200 GeV. The 
dots and solid line shows the central value of the SM and model I prediction, respectively. The 
region between two dashed lines shows the SM prediction: B(B — * K*~^) = (3.25 43) x 10 -5 . 
The shaded band shows the data within 2a errors: B(B —* K*^^) exp = (4.18 ± 0.64) x 10~ 5 . 
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FIG. 2: The Mh dependence of the branching ratio B(B — > K*j) in model II for tan/3 = 4. 
The dot-dashed and solid curve shows the central value of the NLO model II prediction for 
B(B — > K*°j) and B{B — > K*~^), respectively. The region between two dashed lines shows the 
SM prediction: B(B K*~i) = (3.25±^) x 10" 5 . The shaded band shows the same data as 
in Fig[3 
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FIG. 3: The tan/3 dependence of the branching ratio B{B — > K*j) in model II for Mh = 300 
GeV. The curves and bands have the same meaning as in Fig|2 




FIG. 4: The same as FigH but for F K * = 0.38 ± 0.06 instead of F K * = 0.25 ± 0.06. 
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FIG. 5: The Mg dependence of the branching ratio B(B — > K*y) in model III-A. The dot-dashed 
and solid curve shows the central value of the NLO model III-A prediction for B(B — > K 7) 
and B(B — > K*~j), respectively. The region between two dashed lines shows the SM prediction: 
B(B -> K*~i) = (3.25l^) x 10 -5 . The shaded band shows the measured B(B -> K*~i) 
within 2o" errors. 




FIG. 6: The same as FigEJ but for model III-B, i.e. (A ti , X bb ) = (0.5, 22). 
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FIG. 7: The [i dependence of the isospin symmetry breaking Ao-(-KT*7) in the SM and the 
general 2HDM's. The error bar shows the data. For details see the text. 
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FIG. 8: The Mjj dependence of the isospin symmetry breaking Ao-(-K"*7) in the SM and the 
general 2HDM's. The error bar shows the data. For details see the text. 
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FIG. 9: Plots of the M# dependence of CV iSM (mb) (horizontal dots line), AC7 jJVP (dashed curve), 
CV( m &) (solid curve), W 2 \R(K 7) | (short-dashed curve) and W 2 \R(K* 7)! (dot-dashed curve). 
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FIG. 10: Plots of the Mh dependence of the branching ratio B — > /j°7 in the SM (shaded 
band), the model I and II (dashed and dot-dashed curves), and the model III- A (short-dashed 
curve) and III-B (solid curve). The horizontal dots line shows the experimental upper bound 
(at 9%C.L. Q): B{B -» p°j) < 1.2 x KT 6 
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FIG. 11: The same as Fig^Jbut for the decay B — > p 7. 
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FIG. 12: Plots of the angle 7 dependence of the CP asymmetries for B — > p°7 (dashed curve) 
and -B 1 * 1 — > p ± 7 (solid curve) decays in the SM. 
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FIG. 13: Plots of the angle 7 dependence of the CP asymmetry for — > / o ± 7 decay in the SM 
and model III-B for /j, = mi,/2 (dots curves), nib (solid curves ) and 2nib (dashed curves). 
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FIG. 14: The isospin breaking A(pj) vs the CKM angle 7 in the SM and general 2HDM's for 
tan p = 4 and M H = 250 GeV. 
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